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Abstract 

By considering an ensemble of Brownian particles suspended in a heat bath as 
a thermodynamic system with an internal degree of freedom it is possible to ob- 
tain the Fokker-Planck equation for Brownian motion in a temperature gradient, 
by applying the scheme of non-equilibrium thermodynamics. We recover explicitely 
the equations derived in particular by Zubarev and Bashkirov using statistical me- 
chanical and kinetic methods. In addition when the temperature gradient does not 
have an externally imposed magnitude we obtain the differential equation for the 
temperature field, which is coupled to the Fokker-Planck equation. 
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1 Introduction 



In 1965, Nicolis ]IJ considered the problem of Brownian motion in a temperature gradient 
by kinetic methods and derived an expression for the thermal diffusion coefficient in this 
system. A few years later, in 1968, Zubarev and Bashkirov |2) reconsidered this problem 
using the methods developed by Zubarev for treating non-equilibrium systems. They 
explicitely derive the Fokker-Planck equation for this case. In addition they recovered 
Nicolis' expression for the thermal diffusion coefficient. On the other hand Prigogine and 
Mazur ||, || introduced already in 1953 the concept of internal degrees of freedom into the 
formalism of non-equilibrium thermodynamics and pointed out that in this way Fokker- 
Planck like equations could be obtained. Meixner || gave a more general framework for 
the derivation of kinetic equations within this context. 

In this paper we apply the non-equilibrium thermodynamic method of internal degrees 
of freedom to Brownian motion in a temperature gradient and thus derive in a simpler way 
the Fokker-Planck equation for this case. We do not consider the temperature gradient 
as given but establish also the coupled differential equation for the temperature field. In 
section 2 we consider the Brownian particles suspended in a nonuniform heat bath as a 
system with an internal degree of freedom, the velocity u of the particles. We state the 
corresponding Gibbs equation for the total differential of entropy per unit volume and 
give the conservation laws for mass and energy. The entropy balance equation is estab- 
lished in section 3, where we also give expressions for the entropy flux and the entropy 
source strength. In section 4 we follow the principles of non equilibrium thermodynam- 
ics to formulate the linear phenomenological laws for the dissipative fluxes occurring in 
the entropy source strength, also for those pertaining to the internal coordinate space. 
By substitution we then obtain the Fokker-Planck equation and the coupled differential 
equation for the temperature field. Finally, in section 5 we discuss the thermal-diffusion 
regime in which the Fokker-Planck equation may be replaced by a simpler differential 
equation for the density of the Brownian particles. 



2 The Brownian gas as a thermodynamic system with 
an internal degree of freedom 

Consider a dilute 'gas' of Brownian particles of mass m suspended in a heat bath at rest 
with constant mass density pu- The Brownian gas has an internal degree of freedom 
-internal coordinate- namely the velocity u of a Brownian particle. The mass density of 
the Brownian particles of velocity u at position r and time t in the heat bath is p(u, r, t) 
while 

f(u,f,t) = p(u,r,t)/m (2.1) 



2 



is the probability density for the Brownian particles in velocity-coordinate space, which 
may be viewed as an internal thermodynamic variable. 

The system, Brownian gas and heat bath, has mass density p(r, t) 

P = Ph + Pb = Ph + m J f(u,r,t)du, (2.2) 

and energy e(r, t) and entropy s(f,t) per unit of mass. 

Since pu is constant, the Gibbs equation for this system can be written as 

J ffl r 

Sps = -5pe - — J p(u, r, t)5f(u, f, t)du (2.3) 

where 5 denotes the total differential of a quantity and where p(u,f,t) is the chemical 
potential per unit mass of the Brownian gas component with internal coordinate u; T(f) 
is the temperature of the heat bath at position r. 

The chemical potential p(u,r,t) is related to the energy e and entropy s by 

J p(u, f, t)p(u, f, t)du + PhPh = pe- Tps + p (2.4) 

Here pu is the chemical potential of the heat bath and p the system's hydrostatic pressure. 

For the quantities p(u,f,t) = f(u,f,t)/m and e, conservation laws may be written 
down. Let us first deal with conservation of mass. For the rate of change with time of 
f(u, f, t), one may formally write 

d f u df d f ^ df d f 

— = —r ■ — — — ■ Jit = —u ■ — — — ■ Jq (,2.5) 
at or ou or ou 

where represents a flux in velocity space which results from the interaction with the 
heat bath (there is no external field of force acting on the Brownian particles). 

Integrating eq. ( |2.5|) over velocity u we obtain the law of conservation of mass for the 
Brownian particles 

^ = -V • PB v B (2.6) 
where p B is the mass density of the Brownian particles, c.f. eq. ( |2.2|) , and VB(r,t), 

v B (r, t) = p B x J p(u, f, t)udu, (2.7) 
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their mean velocity. 

The law of conservation of energy can be formulated as follows 



where J q represent a heat flux defined in the reference frame in which the heat bath is at 
rest. 



In the next section we shall derive, using the Gibbs equation ( |2.3j ) and the conservation 
laws (|2.5|) and ( |2.8|) , the entropy balance equation for the Brownian gas. 



3 Entropy balance and entropy source strength 

Before deriving an expression for the entropy production in the system described in section 
2, we note that the Brownian gas, i. e. the suspension of Brownian particles in the heat 
bath, may be considered to be a many component ideal solution of 'components' with 
mass density p(u,r,t) and that the thermodynamic potential p(u,r,t) is therefore of the 
form (with k Boltzmann's constant) 

kT 

p(u, r, t) = — lnf(u, r, t) + C(u, f, t) (3.1) 
m 

The potential function C (u, f, t) which depends on the internal coordinate u, can in prin- 
ciple still be a function of the local thermodynamic state variables T(r) and pB{r,t). 

Consider now the gas to be in local equilibrium, i.e. in the state of internal equilibrium 
which would be reached if the system were locally insulated at temperature T(r) and mass 
density ps(r, t). The chemical potential p(u,f,t) and the distribution function f(u,r,t) 
must then satisfy the following requirements: 

1. Since entropy at constant energy and constant density of Brownian particles has in 
that state a maximum as a function of the internal states p(u, r, t), p Leq '(u, r, t) is uniform 
in velocity space and equal to psif^t) the thermodynamic potential of a Brownian ideal 
gas at temperature T(r) and density pB(r,t) 

p Leq iu, f,t)=p B = — [ In^- - -In (3.2) 

m \ m 2 m J 

2. The distribution function f(u) is Maxwellian at local equilibrium and is given by 



4 



f Leq iu,r,t) = exp(ji B - ^ m/kT (3.3) 



If we then substitute eqs. ( |3.2| ) and fl3.3|), for local equilibrium into eq. ( |3.1| ), we find for 

C{u) 

C(u) = l -u> (3.4) 
and therefore that the chemical potential may be written as 

kT 1 

fi(u, f, t) = — lnf(u, f, t) + -u 2 , (3.5) 
m 2 

or equivalently 

kT / \ 

n(u,r,t)= f i B + —ln(f/f Le «) (3.6) 

Let us now calculate the rate of change of entropy per unit volume by differentiating 
eq. (|2.3|) with respect to time 

dps 1 dpe r p(u, f, t) df(u, f, t) 

-dr = T^r- m J— m- (3 - 7) 



Using the conservation laws ( |2.5| ) and ( |2.8| ) as well as eqs. (|3.5| ) and ( p.6|) , expressing 
/j,(u, f, t) in terms of p(u, r, t) = mf(u, r, t), this equation may be written in the form of 
a balance equation 

^ = -VJ s + a, a>0 (3.8) 

where the entropy flux J s and the entropy source strength a, which is a positive quantity 
in accordance with the second law of thermodynamics, are given by 



J s = J'JT -k J f(u, f, t)(lnf(u, f, t) - l)udu (3.9) 
a = - j;T- 2 VT -kjjf—ln (f/f Leq -) du (3.10) 



The modified heat flux J' q 
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J' q = J q — m J -u 2 uf(u,r,t)du (3-H) 

contains a contribution to transfer of heat which is due to the motion of the Brownian 
particles. In deriving eq. (|3.9|) for the entropy flux use has been made of the identity 



f _ df _ d f _ _ r _ df 
k u ■ -—Infdu = k— ■ / uf Infdu — k u- —du 
J or or J J or 



k—- J uf{lnf -l)du (3.12) 



Furthermore, to obtain the second term in eq. (|3.10|) , a partial integration over velocity 



space has been performed, using the fact that the flux vanishes as u — > ±00. 

The entropy source strength a consists of two contributions: the first represents the 
entropy created by heat conduction in the heat bath in the presence of Brownian particles, 
the second, caused by the motion of the Brownian particles in the heat bath, arises so to 
say from diffusion in velocity space, the space of the internal degree of freedom. 



4 Phenomenological relations and the Fokker-Planck 
equation for Brownian motion in a temperature 
gradient 

Following the principles of nonequilibrium thermodynamics 0], the 'linear' phenomeno- 
logical relations -linear relations between the fluxes and thermodynamic forces occuring 
in the entropy production ( |3.10|) - are, since the system is isotropic and assuming locality 



0], 0] in velocity space, 

J' q = -L TT WT/T 2 - J kL TU —ln (f / f Leq ) du (4.1) 

= -L^VT/T 2 - kL u ~ln {f / f Leq ) (4.2) 
where the phenomenological coefficients obey the Onsager-Casimir symmetry relations 

Ut = -L T u (4.3) 
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The assumption of locality in velocity space characterizes the specific physical nature of 
the system considered and is appropriate for Brownian motion. 



With A = Ltt/T 2 a heat conductivity coefficient and defining the friction coefficients 
7 and (3 

7 = Lzr/fT , (3 ee mL m /fT (4.4) 
relations ( [4.1| ) and ( |4.2| ), taking into account relation ( |4.3[ ), become 

4 = _AVT + m / 7 (/ 3+ f|)« (4.5) 

^/vr/r-^+fl) (4.6) 

Considering 7 and /3 to be independent of u in first approximation and substituting 
eq. ([4.6|) into the conservation law ( |2.5| ), we obtain the Fokker-Planck equation for the 
Brownian motion in a heat bath with a non-uniform temperature distribution 

^---u — + (3— (fu+—^-)+^— f— (4 7) 

dt dr du \ m du J T du dr 



This equation is coupled to the differential equation obtained by substituting eq. ( |475l 
into the energy conservation law 



dpe d 2 T d r / _ kT df\ „ , 9 2 T 9 ^ 

"aT = A 1TT ~ ' / / M + — ^ ) du = A 7TT ~ • Pbv b (4.8) 

ar ar z or J \ m ou J or z or 

In writing down eq. ( |4.8| ) we have neglected the small contribution to the heat flux arising 
from the kinetic energy of the Brownian particles. Equation ( |4.7|) coincides, using their 
definitions of constants and variables, with the Fokker-Planck equation derived by Zubarev 
and Bashkirov by statistical mechanical methods, except for one term. This term which 
these authors show to be proportional to the effective volume of the Brownian particle and 
as a consequence of negligible magnitude, can of course, due to its extreme microscopic 
nature, not be found within the framework of a thermodynamic theory. The two equations 
(|4.7| ) and ( |4.8j ) completely specify the coupled evolutions of the temperature field and 
the velocity-coordinate probability distribution of the Brownian particles. However after 
times much larger than the characteristic time the system enters the diffusion and 
thermal diffusion regime for which the evolution is governed by a simpler set of equations. 
We shall discuss this regime in the next section. 
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5 The thermal diffusion regime 



Before discussing the equations which describe the long time behaviour of the Brownian 
particles and the heat bath, we shall derive the equation expressing, for the Brownian 
particles alone, the law of conservation of momentum, an equation which we did not need 
before. This equation will enable us to simplify the equation of motion of the Brownian 
gas for t » (see also in connexion with the developments in this section ref. [|J ch 
X, § 8). 

Using the definition ( |2.7| ) as well as the continuity equation Q2.5| ) for the distribution 
function f(u,r,t) one obtains the following equation of motion for the mean velocity 
v B (f,t): 

PB l[t = ~ V ■P B (r,t)+mJ J^du (5.1) 
Here the hydrodynamic time derivation d/dt for the Brownian gas is defined as 

while its pressure tensor P B is given by 

P B (r,t)=m J f{u-v B )(u-v B )du (5.3) 

Substituting into eq. ( |5.1|) the phenomenological equation (|4.6| ) (which gave rise to the 
Fokker-Planck equation ( [4.7] )) the equation of motion becomes 

^ + P B l V ■ P B (r, t) + 7 VT/T = -(3v B . (5.4) 

As is well-known, in the diffusion regime, i.e. for times long compared to the charac- 
teristic time the Brownian gas, due to collisions with the molecules of the heat bath, 
will reach a state of internal equilibrium (i. eq.). In this state the distribution function 
is approximately a Maxwellian corresponding to a density p B (f,t) and temperature T(r) 
and a non-vamshing mean velocity v B (f,t). It is then given by 

f{u, r, t) ~ f Leq -{u, f, *) = exp{m[p B ~\{u- v B f]/kT} (5.5) 

The pressure tensor P B , eq. ( p.3[ ) then reduces to the gas pressure p B 
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— * — * 

Pb=PbU , p B = p B kT/m, (5.6) 

— * 

with U the unit tensor. At the same time, the inertia term at the left hand side of eq. 

(|5.4j) becomes negligibly small so that this equation, using also eq. (|5.6|) , can be written, 
with Jd = PbVb, as 

J D = -DVp B - D T VT/T (5.7) 

where the diffusion coefficient D and the thermal diffusion coefficient Dt are defined 
respectively as 

D = — R (5.8) 

D T = p B D (l + (5.9) 

The form ( |5.9| ) for Dt agrees with the one obtained by Zubarev and Bashkirov, and earlier 
by Nicolis, by other methods. 



Let us finally discuss the form the entropy production a, eq. ( |3.10| ), takes in the 



diffusion and thermal diffusion regime. We have stated that in this regime f(u,r,t) is 
approximately given by eq. ( |5.5| ) and that in eq. (|5.4|), or the equivalent eq. ([5.1|), the 
inertia term may be neglected, while the pressure tensor becomes diagonal. The following 
relation therefore holds in the diffusion regime 

m J Jn du = VpB (5.10) 
and one then shows with eqs. fl3.3p and ( p.5[ ) that o reduces to 

a = -J q ■ VT/T 2 - J D ■ V PB /pT (5.11) 
Using also eq. (|5.6|) , eq. (|5.11|) may be written as 

a = -I ■ VT/T- - J D ■ {k/m)VpB (5.12) 

PB 

where a new modified heat flux J q has been defined by 
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Jq = Jq+ VbPbJd (5.13) 



Equation (|5.12|) , as does eq. ( |5.11|) , enables one to identify coupled thermodynamic fluxes 



and forces in ordinary space for which Onsager symmetry relations must hold. Let us 
then, with our previous results ( |4.5| ) and ( |5.7| ) for J q and Jd write down the resulting 

phenomenological law for J q . This law is 

J = -AVT - D T T-^ (5.14) 
m p B 

Here the heat conductivity coefficient A at uniform density of the Brownian gas is given 
by 

k n 2 

X = X + -^- (5.15) 
m Dp B 



We note that the coefficients occuring in eqs. ( |5.7| ) and ( |5.14| ) have as expected the correct 
Onsager symmetry. 



ACKNOWLEDGMENTS 

This work has been partially supported by the CICYT of the Spanish Government under 
Grant PB92-0859. 



10 



References 

[1] G. Nicolis, J. Chem. Phys. 43, (1965) 1110. 

[2] D. N. Zubarev and A. G. Bashkirov, Physica 39, (1968) 334. 

[3] I. Prigogine and P. Mazur, Physica XIX, (1953) 241. 

[4] S. R. de Groot and P. Mazur, "Non-Equilibrium Thermodynamics", (Dover Publish- 
ing Co., New York, 1984). See in particular ch. IX, § 8 and X, § 6. 

[5] J. Meixner, Zeitschr. Physik 149, (1957) 624. 



11 



